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Abstract 

Motivated by the first experimental evidence of meson oscillations in the D system, 
we study D° — D° mixing in the Littlest Higgs model with T-parity, we investigate 
its role in constraining the model parameters and its impact on the most interesting 
flavour observables. We find that the experimental data are potentially strongly 
constraining but at present limited by large theoretical uncertainties in the long- 
distance Standard Model contribution to D° — D° mixing. 



Note added 



An additional contribution to the Z penguin in the Littlest Higgs model with T-parity 
has been pointed out in [HE], which has been overlooked in the present analysis. This 
contribution leads to the cancellation of the left-over quadratic divergence in the cal- 
culation of some rare decay amplitudes. Instead of presenting separate errata to the 
present work and our papers [3H6] partially affected by this omission, we have presented 
a corrected and updated analysis of flavour changing neutral current processes in the 
Littlest Higgs model with T-parity in [7]. 



The phenomenon of meson-antimeson oscillation is very sensitive to heavy degrees of 
freedom propagating in the mixing amplitude and, therefore, represents one of the most 
powerful probes of New Physics (NP). In K and B^ s systems the comparison of observed 
meson mixing with the Standard Model (SM) prediction has achieved a good accuracy 
and plays a fundamental role in constraining not only the unitarity triangle but also 
possible extensions of the SM. The evidence for flavour oscillation of the charmed meson 
D°, instead, has been reported only very recently by BaBar [8] and Belle [9], indepen- 
dently. These experimental results have been combined in [10J to which we refer for 
details. Here we just mention that the analysis in \X0\ allows for CP-violation in mixing 
but not in the decay amplitudes where the SM tree-level contributions are expected to 
dominate. CP-violation in D° — D° mixing is also strongly suppressed, in the SM, by 
the small combination of CKM matrix elements V cb V* b . In the presence of NP, however, 
new CP-violating contributions can occur and spoil this feature. 

Combining the new BaBar and Belle measurements of D° — D° mixing parameters 
yields, in particular, an improvement of almost an order of magnitude on the AMp 
constraint, which now reads [10] 

AMd = (11.7 ± 6.8) • 10~ 3 ps^ 1 . (1) 

This first evidence of D° — D° mixing is certainly welcome as, involving mesons with 
up-type quarks, it is complementary to mixing in K and Bd, s systems in providing 
information on NP. In order to discover a signal for NP in AMp, however, one would 
need high confidence that the SM predictions lie well below the present experimental 
limit. 

Unfortunately, the SM calculation of AMp is plagued by long-distance contributions, 
responsible for very large theoretical uncertainties. In fact, unlike B® s — B® s mixing that 
is completely dominated by short- distance effects generated by the top quark, in AMp 
the non-perturbative physics associated with long-distance effects (e. g. propagation of 
light intermediate states) is potentially large and may even dominate over the short- 
distance ones [TT] . 

The short-distance contribution in AM D [T2|[T3]. indeed, is highly suppressed both by 
a factor (m? s — m'^)/M^ v generated by the GIM mechanism and by a further factor (m? s — 
m d)/ m c due to the fact that the external momentum, of the order of m c , is communicated 
to the internal light quarks in box-diagrams. These factors explain why the box-diagrams 
are so small for D mesons relative to K and B^ s mesons where the GIM mechanism enters 
as m 2 c /M^ v and m^/M^ and external momenta can be neglected. Moreover, a recent 
study of AMd has found that NLO (QCD) corrections to short-distance contributions 
interfere destructively with the LO ones, with the net effect (AMp)^ ~ 2-10~ 6 ps~ 1 [14|. 
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Within the SM, then, the short- distance contribution is negligible and a reliable 
theoretical prediction requires an accurate knowledge of the long-distance ones, whose 
estimates follow at present two approaches. The "inclusive" approach, based on the 
operator product expansion (OPE), relies on local quark-hadron duality and on Aqcd/^c 
being small enough to allow a truncation of the series after the first few terms. The charm 
mass, however, may not be large enough for such an approximation. The "exclusive" 
approach, on the other hand, sums over intermediate hadronic states, which can be 
modeled or fit to experimental data. These exclusive contributions, however, need to 
be known with high precision due to cancellations between states within a given SU (3) 
multiplet and, furthermore, the D° is not light enough that its decays are dominated 
by few final states. As a consequence, in the absence of sufficiently precise data, some 
assumptions are required and yield quite model-dependent results. 

While most studies of long-distance contributions find (AMp)^ ^ 10 _3 ps _1 , values 
as high as 10 _2 ps _1 cannot be excluded [ToT - TP?] . The latter estimates, being of the 
order of magnitude of the experimental constraint in ([1]), presently prevent revealing an 
unambiguous sign of NP in D° — D° mixing. In spite of that, in view of future theoretical 
improvements as well as better experimental accuracies, it is certainly interesting to study 
possible NP contributions to AMp in specific extensions of the SM. It is important to 
note that NP contributions appear in box-diagrams with internal new heavy particles 
and, therefore, are of short-distance only. In predictive NP models, a reliable calculation 
is then possible and its remaining uncertainty is dominated by the parameters of the 
model itself. 

The aim of the present Letter is to study the phenomenon of D° — D° mixing in the 
Littlest Higgs model with T-parity (LHT) and to investigate its impact on our previous 
LHT flavour analyses [3|fT8]. 

The LHT model [191120] belongs to the class of the so-called Little Higgs models 
(LH) [21], whose basic idea for solving the little hierarchy problem is the interpretation 
of the Higgs as a pseudo-Goldstone boson of a spontaneously broken global symmetry. 
Diagrammatically, the quadratic divergences that affect the Higgs mass, within the SM, 
are canceled by the contributions of new heavy particles having the same spin-statistics as 
the SM ones and masses around 1 TeV. In the LHT model, a discrete symmetry called T- 
parity is added, in order to satisfy the electroweak precision constraints [22], by avoiding 
tree-level contributions of the new heavy gauge bosons and restoring the custodial SU(2) 
symmetry. Under T-parity particle fields are T-even or T-odd. The T-even sector 
consists of the SM particles and a heavy top T+, while the T-odd sector contains heavy 
gauge bosons (W^, Z H , A H ), a scalar triplet (<&) and the so-called mirror fermions, i.e. 
fermions corresponding to the SM ones but with opposite T-parity and 0(1 TeV) mass. 
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Mirror fermions are characterized by new flavour violating interactions with SM fermions 
and heavy gauge bosons, which involve two new unitary mixing matrices in the quark 
sector, analogous to the Cabibbo-Kobayashi-Maskawa (CKM) matrix Vckm- They are 
Vud and Vh u , when the SM quark is of down- or up- type respectively, and they satisfy 
VhuVhcI = ^ckm [23] ■ A similar structure is valid for the lepton sector as discussed 
in details in jl]. It is important to recall two important features of the LHT model 
in order to understand its role in flavour physics. The first is that, because of these 
new mixing matrices, the LHT model does not belong to the Minimal Flavour Violation 
(MFV) class of models whether constrained |24j or general [25J and significant effects 
in flavour observables are possible. The second is that no new operators, and no new 
non-perturbative uncertainties, in addition to the SM ones appear in the LHT model. 

Extensive flavour physics analyses in the LHT model have been recently performed 
in both the quark [31 El HSl E3J ESI [27] and lepton sector [SEE]. In particular, D° - D° 
mixing has been studied in [Hl[23], before it was experimentally observed. Motivated 
by the improved experimental constraint on AMd [BHTU] we update and extend here the 
LHT analysis of D° — D° mixing. 

As discussed above, at present the large SM long-distance uncertainties prevent to 
reveal an unambiguous NP contribution to AMd- We choose, therefore, to disentangle 
our analysis from the large SM uncertainties. To this end, we consider only the LHT 
contribution to AMd and determine the range of values that it can assume once the 
known flavour constraints are imposed as in [21 [T8] . Once the SM uncertainties are 
significantly reduced, our strategy can be pushed further to use the experimental AMd 
measurement to constrain the parameters of the LHT model. Moreover, if the smaller 
SM upper bounds (AMd)sm ^ lCT 3 ps -1 are confirmed, it will be legitimate to neglect 
the SM contributions and to use the experimental AMd measurement as a constraint 
for the LHT contribution only. 

Meson-antimeson mixing in the LHT model is discussed in details in JT8JI23] where 
the separate contributions of T-even and T-odd sector to the off-diagonal dispersive 
matrix elements M{ 2 {i — K,d,s for K and B^s systems) are explicitely given. The D 
system presents a difference since it involves external SM up-type quarks and therefore 
the T-even T + cannot run in the loop, so that the T-even contribution reduces to the 
SM one. In the T-odd sector, both down-type mirror quarks, together with the charged 
gauge bosons Wjj, and up-type mirror quarks, together with the neutral gauge bosons 
Zh,Ah, contribute. 

Due to the near equality of up- and down-type mirror quark masses the formula 
for the AC = 2 effective Hamiltonian can straightforwardly be obtained from the one 
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describing AS = 2 transitions calculated in [Hl[23], yielding 



[n eS (AC = 2)] odd = ^M^oJ^^^Mzi, h) {uc) v . A {uc) v . A . (2) 



,2 //\#2 



Here, Zi = m Hi / M^ Vh with mm denoting the mass of the z-th mirror quark doublet, and 
the function Fh has been determined in jT8.,23j. The only difference with respect to 



[H e e(AS = 2)] dd is the fact that now the relevant quark mixing is given by the Vhu 
matrix, leading to the combination ^• D ' ) \_ 

tiP) \r*iu\ric ( n\ 

St — V Hu V Hu ■ \ 6 ) 

The QCD correction t]d can be approximated by [29J 

Vd ~ r] 2 = 0.57 ± 0.01 . (4) 

The mirror quark contribution to the off-diagonal element M® 2 of the neutral D-meson 
mass matrix is then found to be 

(M3 o4d ^|(Mg) od Je-^ 



^FlB D m D0 M 2 w V - V2 d D) *ef r F H {z u z 3 ) . (5) 

Our definition of the phase <f)£> follows from 

(Mg)* = (D°\n cS {AC = 2)\D°) = | (Mg) | e 2ife , (6) 

where we note that it is (Mj^)* and not M^, as sometimes found in the literature, 
that appears on the 1. h. s. of ([6]). We stress that the theoretical uncertainty on the 
LHT contribution, being of short- distance origin only, comes from the non-perturbative 
uncertainties in the decay constant Fd and the .B-parameter Bp, in addition to the new 
LHT parameters scanned over in the analysis. For Fd we use the recent experimental 
determination by the CLEO-c collaboration [30J that turned out to be in agreement 
with recent lattice calculations [3T1 [32] and of comparable precision. Concerning the 
parameter B D , we consider the result of the most recent (quenched) lattice calculation 
[31], compatible within quite large uncertainties with an older lattice determination [33J. 
Their numerical values are collected together with the other inputs in Table [U where 
some numbers have been updated with respect to [3lfl8]. 

Having at hand all the LHT formulae for meson oscillations in K, D and B^ s systems 
and rare K and B decays presented in [3|[T8] and here, we will now investigate the impact 



1 We would like to caution the reader that in [18] the indices of the Vh u elements have been erroneously 
interchanged. 
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Gj? = 1.16637 • 10 GeV 


-2 


AM K = 3.483(6) • 10" 15 GeV 


Mw — 80 425(38) GeV 




AM C 


= 0.508(4)/ps 


[35] 


a = 1/127.9 
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m = 


1.32(32) 


rrn 

m] 


mjfo = 497.65(2) MeV 




m = 


0.47(5) 


[42] 


m D o = 1.8645(4) Lxev 




m = 


0.57(1) 




m Bd = 5.2794(5) GeV 




VB = 


0.55(1) 


[29] 


m Bs = 5.370(2) GeV 




m c = 


= 1.30(5) GeV 




\e K \ = 2.284(14) • 10~ 3 


[33] 


m t = 


= 161.7(20) GeV 




S^ Ks = 0.675(26) 


[35] 









Table 1: Values of the experimental and theoretical quantities used as input parameters. 



of the measurement ([I]) and of the constraints on \M^\ and (ftp derived in [10] on the 
parameters of the LHT model and our results presented in 0,[18]. To this end we will 
consider two frameworks for the SM contributions to D° — D° mixing: 

Framework X 

The SM contribution to M® 2 is set to zero so that the constraint on \M® 2 \ and the 
phase 4>d shown in the lower left plot in Fig. 2 of [10] is directly applied to the LHT 
contribution. 

Framework Y 



The SM contribution is allowed to vary within its large uncertainties as done in [10J and 
the general constraint on the NP contribution shown in the lower right plot in Fig. 2 
of [10] is applied to the LHT model. 

In Fig. [I] we show the predictions of the LHT model for \M^\ and 20d obtained in 
a general scan (blue points) over the parameters of the model in comparison with the 
allowed la ranges 2 ) derived in [TP] . If we allowed for the 2a ranges instead, there would 
be almost no restrictions on the Little Higgs parameter space from D° — D° mixing, 
i. e. there would be no visible difference between the plots with the points allowed in 
frameworks X and Y. We therefore restrict ourselves to the la ranges where the effect is 



2 In practice the constraints derived in |10) have been implemented in our analysis, in the X and Y 
Frameworks respectively, approximating the la ranges as: X) 0.0025 ps -1 < < 0.0125 ps^ 1 and 

2|<M < 50°; Y) \Mg\ < 0.005 ps^ 1 or (2|^| < 25° and \Mg\ < 0.02 ps" 1 ). 
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Figure 1: \M^ 2 \ versus 2<p£, from a general scan over the LHT parameters, compared to 
the probability density function derived in [10], for the Framework X (left) and Y (right). 
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Figure 2: The same as in Fig. U\ but for two LHT specific scenarios: K -Scenario (brown 
points) and B s -Scenario (green points). 



quite pronounced. If in the future the 2a ranges come down to where there are now the 
la ranges, D° — D° mixing and in particular its CP-violating phase will put significant 
restrictions on the Little Higgs parameter space. Fig. [2] shows analogous results in two 
specific parameter scenarios identified in [3],[T8]: the iT-scenario (brown points) and B s - 
scenario (green points) that lead to large departures from the SM in K and B decays, 
respectively. Finally in Figs. IMSl we show the impact of the experimental D° — D° 
constraint on the most interesting results found in [5l[T5] . 

From the inspection of Figs. [Q15] and the comparison with our previous results [3|[T8] 
we learn that: 

• The D° — D constraint is much weaker in the Framework Y, due to very large 
long-distance uncertainties present in the SM. In the Framework X the latter are 
only present in FdV Ad and, as seen in Figs. [T]l2l the impact of the D° — D° 
constraint on the points satisfying all remaining observables is rather significant. 

• We observe from Fig. [2] that whereas the ii'-scenario is practically excluded by the 
D° — D° mixing data in the Framework X, the impact on the S s -scenario is only 
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Figure 3: Bt(Kl — > tt uu) as a function of Br(K + — > 7i + uu), after applying the la- 
constraint on D° — D° mixing within the Framework X (left) and Y (right). The shaded 
area represents the experimental la-range for Br(K + — > 7i + uu). The GN-bound U^j is 
displayed by the dotted line, while the solid line separates the two areas where BriK^ — > 
^vv) is larger or smaller than Br(K + — > -r^vv). 



Br (K L ^7r°e*e~ / jx* jjT ) Br (K L ^7T°e*e~ / u' W ) 




Figure 4: BriK^ — > 7r°e + e~) (upper curve) and BriK^ — > 7i°fi + fi~) (lower curve) as 
functions of Br(Ki — > 7i uu), after applying the la-constraint on D° — D° mixing within 
the Framework X (left) and Y (right). 
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Figure 5: Br(Ki — > 7t°i>b>) as a function of S^<j>, after applying the la-constraint on 
D° — D° mixing within the Framework X (left) and Y (right). 
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moderate in both SM frameworks. Therefore, the impact of the D° — D° constraint 
turns out to be significantly larger on K decays than B decays. The reason is that 
both K decays and D° — D° mixing describe transitions between the first two 
quark generations, thus involving the same combinations of elements of Vnd and 
Vhui respectively. Now, as Vnd and Vhu are related via V Hu = V Hd V<l KM and 
Vckm — 1, it approximately turns out that Vh u — Vna,- Therefore, the observed 
correlation between K and D physics in indeed expected within the LHT model. 

• As shown in Figs. SO in the case of the Framework Y and a general scan over 
the LHT parameters, very large departures from the SM expectations for rare K 
decays and are possible. These plots, in fact, are qualitatively similar to those 
presented in [31 ITS]. 

• On the other hand, if mirror fermion contributions describe the full D° — D° mixing 
as supposed in the Framework X, the enhancements of rare K decay branching 
ratios are significantly smaller than those found in [3], although they can still be 
substantial. For instance Bt(Kl —> 7t°i>u) can be larger by a factor 5 relative to 
the SM prediction. On the other hand the CP-conserving decay Br(K + — > tt + vv) 
is less affected by the D° — D° mixing constraint. 

• Finally we observe that a large phase (pr> can be generated, signaling the possibil- 
ity of sizeable CP-violating effects in the D meson system within the LHT model. 
Quantitative predictions for CP-violating observables in the D system would how- 
ever require a much more detailed analysis which is beyond the scope of the present 
Letter. 

The main message of our paper is that the present data on D° — D° mixing put 
already significant constraints on the predictions of the LHT model for K and B decays. 
However, without a consistent improvement in the estimate of the SM contribution to 
D° — D° mixing, the role of the D system in constraining the parameters of the LHT 
model as well as other extensions of the SM will be limited, even if the accuracy of the 
data improves. The situation is more promising in the case of CP-violation in the D 
meson system, where due to the absence of SM contributions much cleaner predictions in 
a given NP model can be made. On the other hand, useful constraints on the parameter 
space of a given NP model can only be obtained once the data significantly improve. 

Acknowledgements 

We thank Ikaros Bigi for an illuminating and informative discussion on the physics of 
D° — D° oscillations. This research was partially supported by the Cluster of Excel- 



8 



lence 'Origin and Structure of the Universe' and by the German Bundesministerium fur 
Bildung und Forschung under contract 05HT6WOA. 



References 

[1] T. Goto, Y. Okada and Y. Yamamoto, Phys. Lett. B 670 (2009) 378 
|arXiv:0809.4753l [hep-ph]]. 

[2] F. del Aguila, J. I. Illana and M. D. Jenkins, JHEP 0901 (2009) 080 
larXiv:0811.289l1 [hep-ph]]. 

[3] M. Blanke, A. J. Buras, A. Poschenrieder, S. Recksiegel, C. Tarantino, S. Uhlig 
and A. Weiler, JHEP 0701, 066 (2007) |arXiv:hep-ph/061029 8|. 

[4] M. Blanke, A. J. Buras, B. Duling, A. Poschenrieder and C. Tarantino, JHEP 
0705, 013 (2007) [arXiv:hep-ph/0702136] . 

[5] M. Blanke, A. J. Buras, S. Recksiegel, C. Tarantino and S. Uhlig, JHEP 0706, 082 
(2007) [ar Xiv:0704.3329l [hep-ph]]. 

[6] M. Blanke, A. J. Buras, S. Recksiegel and C. Tarantino. larXiv:0805.43 93 [hep-ph]. 

[7] M. Blanke, A. J. Buras, B. Duling, S. Recksiegel and C. Tarantino, [arXlvi0906. 5454 
[hep-ph] . 

[8] B. Aubert et al. [BABAR Collaboration], Phys. Rev. Lett. 98 (2007) 211802 
|a?Xlv":hep-ex/0703020] . 

[9] M. Staric et al. [Belle Collaboration], Phys. Rev. Lett. 98 (2007) 211803 
[arXiv:hep-ex/0703036] . K. Abe et al. [BELLE Collaboration], larXiv: 0704.10001 
[hep-ex] . 

[10] M. Ciuchini, E. Franco, D. Guadagnoli, V. Lubicz, M. Pierini, V. Porretti and 



L. Silvestrini, [arXiv:hep-ph /0703204, 
[11] For detailed discussions and recent reviews see: G. Burdman, 



arXiv:hep-ph/9407378. G. Raz, Phys. Rev. D 66 (2002) 057502 



|arXiv : hep-ph/ 02 05 1 1 3) . S. Bianco, F. L. Fabbri, D. Benson and I. Bigi, 



Riv. Nuovo Cim. 26N7 (2003) 1 arXiv:hep-ex/0309021|. G. Burdman and 



I. Shipsey, Ann. Rev. Nucl. Part. Sci. 53 (2003) 431 | arXiv:hep-ph/0310076| . 
A. A. Petrov, Int. J. Mod. Phys. A 21, 568 (2006) |arXiv:hep-ph/0611361] . Y. Nir, 



arXiv:hep-ph/0703235 



9 



[12] H. Y. Cheng, Phys. Rev. D 26 (1982) 143. 



[13] A. Datta and D. Kumbhakar, Z. Phys. C 27 (1985) 515. 

[14] E. Golowich and A. A. Petrov, Phys. Lett. B 625 (2005) 53 
[arXiv:hep-ph/0506185] . 

[15] I. I. Y. Bigi and N. G. Uraltsev, Nucl. Phys. B 592 (2001) 92 
[arXiv:hep-ph/00050"8~9] . 

[16] J. Donoghue, E. Golowich, B. Holstein and J. Trampetic, Phys. Rev. D 33, (1986) 
179. L. Wolfenstein, Phys. Lett. B 164, (1985) 170. P. Colangelo, G. Nardulli and 
N. Paver, Phys. Lett. B 242, (1990) 71. T.A. Kaeding, Phys. Lett. B 357, (1995) 
151 [arXiv:hep-ph/9505393 | . A. A. Anselm and Y. I. Azimov, Phys. Lett. B 85, 
(1979) 72. 

[17] A. F. Falk, Y. Grossman, Z. Ligeti and A. A. Petrov, Phys. Rev. D 65, (2002) 
054034 larXiv:hep-ph/0110317| . A. F. Falk, Y. Grossman, Z. Ligeti, Y. Nir and 
A. A. Petrov, Phys. Rev. D 69, (2204) 114021 |arXiv:hep-ph/0402204] . 

[18] M. Blanke, A. J. Buras, A. Poschenrieder, C. Tarantino, S. Uhlig and A. Weiler, 
JHEP 0612 (2006) 003 |arXiv:hep-ph/0605214| . 

[19] N. Arkani-Hamed, A. G. Cohen, E. Katz and A. E. Nelson, JHEP 0207 (2002) 
034 |arXiv:hep-ph /0206021 



[20] H. C. Cheng and I. Low, JHEP 0309 (2003) 051 |arXiv:hep-ph703 08199|; 
JHEP 0408 (2004) 061 | arXiv:hep-ph/0405243 |. I. Low, JHEP 0410 (2004) 067 
|arXiv:hep-ph/0409025] . 

[21] N. Arkani-Hamed, A. G. Cohen and H. Georgi, Phys. Lett. B 513 (2001) 232 
[arXiv:hep-ph/0105239] . 

[22] J. Hubisz, P. Meade, A. Noble and M. Perelstein, JHEP 0601 (2006) 135 
[arXiv:hep-ph/05060"42] . 



[23] J. Hubisz, S. J. Lee and G. Paz, JHEP 0606 (2006) 041 |ar Xiv:hep-ph/051 2169|. 

[24] A. J. Buras, P. Gambino, M. Gorbahn, S. Jager and L. Silvestrini, Phys. Lett. 
B 500 (2001) 161 |arXiv:hep-ph/0007085] . A. J. Buras, Acta Phys. Polon. B 34 
(2003) 5615 | arXiv:hep-ph/0310208] . M. Blanke, A. J. Buras, D. Guadagnoli and 
C. Tarantino, JHEP 0610 (2006) 003 [arXiv:hep-ph/0604057] . 



10 



[25] G. D'Ambrosio, G. F. Giudice, G. Isidori and A. Strumia, Nucl. Phys. B 645 



(2002) 155 arX iv:hep-ph/0 207036|. 



[26] M. Blanke, A. J. Buras, A. Poschenrieder, S. Recksiegel, C. Tarantino, S. Uhlig 
and A. Weiler, Phys. Lett. B 646 (2007) 253 |arXiv:hep-pn/0609284]. 



[27] H. Hong-Sheng, Phys. Rev. D 75 (2007) 094010 |arXiv:hep-p h/0703067|. 



[28] S. R. Choudhury, A. S. Cornell, A. Deandrea, N. Gaur and A. Goyal, Phys. Rev. 
D 75 (2007) 055011 |a rXiv:hep-ph/0612327] . 

[29] A. J. Buras, M. Jamin and P. H. Weisz, Nucl. Phys. B 347 (1990) 491. J. Ur- 
ban, F. Krauss, U. Jentschura and G. Soff, Nucl. Phys. B 523 (1998) 40 
|arXiv:hep-ph/9710245] . 

[30] M. Artuso et al. [CLEO Collaboration], Phys. Rev. Lett. 95 (2005) 251801 
[arXiv:hep-ex/0508057] . 

[31] H. W. Lin, S. Ohta, A. Soni and N. Yamada, Phys. Rev. D 74 (2006) 114506 
[arXiv:hep-lat/0607035] . 



[32] C. Aubin et al, Phys. Rev. Lett. 95 (2005) 122002 arXiv:hep-lat/0506030|. 



[33] R. Gupta, T. Bhattacharya and S. R. Sharpe, Phys. Rev. D 55 (1997) 4036 
[arXiv:hep-lat/96110~23~| . 

[34] W. M. Yao et al. [Particle Data Group], J. Phys. G 33 (2006) 1. 

[35] E. Barberio et al. [Heavy Flavor Averaging Group (HFAG) Collaboration], 



arXiv:0704.3575 [hep-ex]; http:/ /www. slac.stanford.edu/xorg/hfag/ 



[36] E. Blucher et a/., arXiv:hep-ph/0512039 



[37] M. Bona et al. [UTfit Collaboration], JHEP 0603 (2006) 080 
[arXiv:hep-ph/0509219] ; Phys. Rev. Lett. 97 (2006) 151803 
|arXiv:hep-ph/0 605213|; http://utfit.romal.infn.it 



[38] A. Abulencia et al. [CDF Collaboration], Phys. Rev. Lett. 97 (2006) 242003 
[arXiv:hep-ex/0609040] . 

[39] V. M. Abazov et al. [DO Collaboration], Phys. Rev. Lett. 97 (2006) 021802 
[arXiv:hep-ex/0603029] . 



[40] S. Hashimoto, Int. J. Mod. Phys. A 20 (2005) 5133 |arXiv:hep-ph/0411126| . 



11 



[41] S. Herrlich and U. Nierste, Nucl. Phys. B 419 (1994) 292 |arXiv:hep-ph/9310311| . 



[42] S. Herrlich and U. Nierste, Phys. Rev. D 52 (1995) 6505 |arXiv:hep-ph/ 9507262 ] ; 
Nucl. Phys. B 476 (1996) 27 |arXiv:hep-ph/9604330] . 

[43] Y. Grossman and Y. Nir, Phys. Lett. B 398 (1997) 163 | arXiv:hep-ph/9701313| . 



12 



